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[j^ I Abstract 

We consider the polarization behaviour of radio waves propagating through 

an ultrarelativistic highly magnetized electron-positron plasma in a pulsar mag- 

netosphere. The rotation of magnetosphere gives rise to the wave mode coupling 

^ I in the polarization-limiting region. The process is shown to cause considerable 

circular polarization in the linearly polarized normal waves. Thus, the circular 
polarization observed for a number of pulsars, despite the linear polarization of 
fvq . the emitted normal waves, can be attributed to the limiting-polarization effect. 
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0\ ■ 1 Introduction 

::a: 

f^ , Pulsar radiation is known to have a high linear polarization. Although the cir- 

cular polarization of the average pulse profiles is typically less than the linear 
one, for the most of pulsars the circular polarization is also significant (for recent 
review, see Han et al. 1998 and references therein). The circular polarization 

5^ . tends to peak near the centre of the profile. In many pulsars the circular po- 

larization reverses its sense near the pulse centre, while in the others the sense 
remains constant throughout the pulse (Radhakrishnan & Rankin 1990). The 
^\f ' observed circular polarization of pulsar radiation is commonly attributed either 

to the propagation effects or to the emission mechanism. 

The magnetic field pattern of a pulsar is generally considered to be that 
of a rotating dipole. The magnetosphere contains an electron-positron plasma 
outflowing along the open magnetic field lines with Lorcntz- factors 7 ~ 100. 
So pulsar radiation, believed to originate in the open field line tube, propa- 
gates through the highly magnetized relativistic plasma. The propagation of 
radio waves in a relativistically moving medium was analyzed by Lee & Lerche 
(1975), Harding & Tademaru (1981) However, they considered rather abstract 
model of a medium, with the spatial dispersion being ignored. A qualitative 
consideration of wave polarization behaviour in application to a more realistic 
pulsar model allowed to attribute a number of characteristic features of the ob- 
served polarization pattern to the propagation effects inside the magnetosphere 
(Cheng & Ruderman 1979, Stinebring 1982, Barnard 1986, von Hoensbroech et 
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al. 1998). 

Given that the plasma is embedded in an infinitely strong magnetic field, 
the normal wave modes are linearly polarized. One of them, namely the ordi- 
nary mode, is polarized within the same plane as the wave vector, k, and the 
magnetic field, b; the other one, called the extraordinary mode, is polarized 
perpendicularly to the this plane. In the emission region, the scale length for 
beats between the wave modes is much less than the scale length for change 
in the plasma parameters. So in the vicinity of the emission origin the normal 
waves propagate independently, their polarization plane being adjusted to the 
local orientation of the magnetic field. Along the trajectory, the angle 6 between 
the wave vector and the magnetic field increases due to magnetic line curvature 
and the plasma density decreases due to open field line tube widening, the scale 
length for beats increasing. Eventually the plasma density decreases enough 
and the medium does not infiuence the waves; their polarization remains now 
fixed. Thus, the emergent polarization is set up by the processes in the so- 
called polarization-limiting region, where the geometrical optics approximation 
is violated (Budden 1952). The polarization-limiting radius, Tp, can be found 
from the relation: 

^An(rp)s(rp) ^ 1, (1.1) 

where An is the difference between the refractive indices of the wave modes 
considered, lu the wave frequency, s the scale length for change in the plasma 
parameters, s ^ rp. In the highly magnetized relativistic plasma of number 
density N and Lorentz-factor 7, the refractive indices for the ordinary and 
extraordinary modes are known to be, respectively (Melrose & Stoneham 1977); 

uj?, sin^ 9 
- -' ^ n, = l, (1.2) 
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with ujp = y jy^^ being the plasma frequency, (3 the plasma velocity in units 
of c. 

In application to pulsars, the polarization- limiting radius was estimated by 
Cheng & Ruderman (1979), Stinebring (1982) and Barnard (1986). Note that 
if a normal wave were emitted in the field line plane of a non-rotating dipolar 
magnetic field, the polarization plane would not vary along the trajectory at all; 
the normal modes then would propagate independently preserving their initial 
polarization states. As the wave propagates in the rotating magnetosphere, the 
wave vector acquires a tilt to the ambient magnetic line planes, so that the kx b- 
plane turns along the trajectory. Another consequence of rotation, namely the 
magnetic line sweepback, also leads to rotation of the k x b-plane, however, this 
effect is less significant. As the scale length for the wave mode beats increases 
sufficiently, the polarization of a normal wave has no time to follow the variations 
of the ambient magnetic field and the wave mode coupling holds. Cheng & 
Ruderman (1979) pointed out that the propagation of the linearly polarized 



normal waves in the region where the geometrical optics fails, results in elliptical 
polarization of the emergent radiation. The effect is likely to account for the 
circular polarization observed in pulsars that exhibit predominantly one sense 
of circular polarization across the pulse (Radhakrishnan & Rankin 1990). 

The influence of pulsar rotation on the observed polarization was also consid- 
ered by Blaskiewicz et al. (1991). These authors investigated the contribution 
of rotation to the velocity of the particles emitting curvature radiation, with 
the propagation effects being neglected. However, it is the wave propagation 
through the magnetosphere that sets up the characteristics of outgoing radia- 
tion. Indeed, whatever the emission mechanism, only the waves corresponding 
to the modes allowed by the magnetospheric plasma can propagate and ulti- 
mately escape from pulsars. 

Provided that the plasma with the various distribution functions for elec- 
trons and positrons is embedded in the finite magnetic field, the normal waves 
propagating at the small angle to the field should be circularly polarized. Pro- 
ceeding from this, Cheng & Ruderman (1979), von Hocnsbroech it et al. (1998) 
proposed to explain the observed circular polarization of pulsar radiation by 
the dispersive properties of the magnetospheric plasma. However, well within 
the pulsar magnetosphere the magnetic field strength is so high that the critical 
angle for the circularly polarized normal modes appears to be too small. 

The present paper deals with a more detailed consideration of wave polariza- 
tion behaviour in the polarization-limiting region. In Sect. 2 we treat the equa- 
tions describing the waves which propagate through an ultrarelativistic highly 
magnetized plasma rotating together with the magnetic field. In Sect. 3 these 
equations are applied to the pulsar magnetosphere. The wave mode coupling in 
the polarization-limiting region is found to produce a significant circular polar- 
ization in the initially linear normal waves. The results are summarized in Sect. 
4. 

2 Basic equations 

Consider an ultrarelativistic electron-positron plasma rotating together with an 
infinitely strong inhomogeneous magnetic field. For the waves of frequency lo 
the wave fields E, B are described by Maxwell's equations: 

_ „ ii-^ „ 47r . 

VxB = E+— ji, 

c c 

VxE=— B, (2.1) 

c 

—iuje(n^ — n^) + divji = 0. 
Here ji is the linearized current density caused by the waves, 

ji =e[vo(n+ -nj;)+no(v+ - Vj^)], (2.2) 



with vq, no being the unperturbed particle velocity and number density, re- 
spectively, Vi , Hi the small perturbations of these quantities for electrons and 
positrons. Although the distribution functions for electrons and positrons are 
generally believed to be slightly different, the difference is insignificant in our 
case. Therefore we assume that Vq and no are equal for both the particle species. 
In the infinitely strong magnetic field, the particle motion can be treated in 
terms of the mechanical bead-on-a-wire model. The particle moves along the 
magnetic field line which rotates at the angular velocity il. The particle velocity 
in the laboratory frame may be written as 

v = rJxr + Wfcb, (2.3) 

where b is the unit vector along the magnetic field, vt the velocity along the 
field line. Substituting Eq. (2.3) into the Lagrangian 



L = -mc^^/1 2" + -(A • v) - eip, 



with A and (p being, respectively, the vector and scalar potentials of the wave, 

in X rl 
one can get the equation of motion. To the first order in - — ^ — - we obtain: 



m-f^^ = ±e 
at 



(n X r) X B 
E H 



(2.4) 



Here the curvature radius of the magnetic field lines was assumed to be much 
higher than the wavelength; indeed, this is a good approximation in our case. 
The equation of motion (2.4) implies that the inertial forces introduced by the 

rotation are small as O ( 2 — I • Note that the left-hand side of Eq. (2.4) 

contains the total derivative, -^ = —iuj + vq • V. 

Equations (2.1), (2.2) and (2.4) yield the self-consistent description of the 
wave fields and the plasma particle motion in these fields. The refractive indices 
for the waves in the polarization-limiting region equal unity to within -j^ <^ 1. 
Therefore the waves propagate almost straight-line and one can choose the three- 
dimensional Cartesian system with the z-axis aligned with the wave vector. 
Then all the perturbed quantities should depend only on the z-coordinate. On 
the basis of Eqs. (2.1) and (2.4) one can write the component equations: 
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We consider the waves propagating along the z-axis in positive direction. 
Since the refractive indices of the waves considered are very close to unity, the 
spatial dependence should be close to exp [i—zj. Then one can present the field 
components in the form: 

Efj^ = a^cxp[i—z), (2.6) 



c 
with the amplitudes a^ varying slowly: 

da, ^^ a^ 
dz c 

Of course, the rest of the perturbed quantities can be presented similarly. The 
scale length for change in the medium parameters, namely in b, Vq and no, also 
exceeds the wavelength essentially. Then the set of equations (2.5) is reduced 
to the form: 



where 



^ = -iR[(b^ + qyfa^ + (b^ + qy)(by - qx)ay], 
-^ = -iR[{bx + qy){by - qx)ax + {by - qxYay], 



b X (ri X r) Wp 



(2.8) 
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with ujp = y ^^^ being the plasma frequency, P^ the z-component of the 
plasma velocity vq in units of c. Above we took into account that the waves 
considered are quasi-transverse ones, 1 — n <C 1, and, correspondingly, the z- 
component of the wave electric field is small, a^ <C Oj;, Oj,. With n given by Eq. 
(1.2), this is surely valid at the condition (1.1). 

Neglecting the rotation and the z-dependence of unperturbed plasma pa- 
rameters we come to the homogeneous problem. Then the solutions are easy to 
be found in the form: 

a,{z) oc exp ( — « — (1 — n)z ) . (2-9) 

Setting the determinant of the system equal to zero immediately yields the 
customary refractive indices (1.2) for the ordinary and extraordinary modes in 
the homogeneous highly magnetized ultrarelativistic plasma. 



3 Polarization transfer in the rotating magneto- 
sphere 

3.1 The main features of wave propagation 

First we examine the set of equations (2.8) to outline the characteristic features 
of wave polarization behaviour along the trajectory in pulsar niagnetosphere. 
As long as the plasma density is sufficiently high, so that Rz ^ 1, geometrical 
optics approximation is valid. Then the wave field amplitudes can be presented 
in the form: 

a^=a|,")exp[G(z)], (3.1) 

with G{z) being as large as Rz. Substituting Eq. (3.1) into Eq. (2.8) and setting 
the determinant of the system equal to zero one can find, to the first order in 
{Rz)-\ 

dG\ 

— j = -iR[{b:, + qyf + {by - q,)\ 

/ e 

with the subscripts o and e referring to the ordinary and extraordinary modes, 
respectively. Note that the terms q^ and qy arc introduced by the magnetosphere 
rotation, so that in the corotating frame g^; = g^ = 0. Comparison of Eqs. (3.1), 
(3.2) with Eqs. (1.2), (2.9) corresponding to the homogeneous plasma then leads 
to a well-known result: the wave vector in a weakly inhomogcneous medium 
equals, in each point of the trajectory, to that in the homogeneous medium with 
the same parameters, so that a^ oc exp [—1^/^(1 — n{z))dz\. With Eq. (3.2) 
we obtain: 

( ax_\ ^ b^ +qy I ax \ _ by - q^ 

UfA by-'i^' V4"7e ^^ + '^^' 
Hence, in the corotating frame the polarization of normal waves follows the slow 
variations of the ambient magnetic field. 

Since the plasma density decreases along the trajectory, N oc z~^, geomet- 
rical optics ultimately fails. In the opposite limit, Rz <C 1, wave propagation 
is no longer influenced by the plasma and, obviously, wave polarization remains 
fixed (cf. Eq. (2.8)). Thus, the emergent polarization is formed in the so-called 
polarization-limiting region, where Rz ~ 1. In this region, the scale length for 
wave mode beats becomes so large that the polarization of a normal wave has 
no time to follow the variation of the ambient magnetic field, so that the wave 
mode coupling holds. 

As follows straightly from Eq. (2.8), if the magnetosphere were non-rotating 
one, the normal waves emitted in the plane containing the field lines would 
propagate independently throughout the trajectory, the polarization plane re- 
maining constant. Now we are interested in the wave mode coupling introduced 



by the slow magnetosphere rotation. The rotation causes the wave vector tilt 
of ^ zJTh to the ambient field line planes along the trajectory; here tl = -r; 
is the light cylinder radius. In addition, the dipolar magnetic field structure 
should be distorted. Indeed, the poloidal current flow, j ~ -^ti^' ^^^^ ^p ^^^^S 
the poloidal magnetic field strength (Goldreich & Julian 1969), gives rise to the 

azimuthal component of the magnetic field, _B^: (rotB)p r^ —j- ^ -^ji where d 

3/2 

is the open field line tube diameter, d — 1I2 • '^^'^ field line twist then can be 



r 



R / \ 3/2 
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estimated as -^ ^ \^] ■ Thus, the sweepback of magnetic lines appears 
to be less significant than the rotation effect. Therefore below we consider the 
rotation of purely dipolar magnetic field. 

Let 5 be the opening half-angle of the emission cone. Consider the wave 
emitted along the z-axis at the angle 5 to the magnetic axis, the latter lying 
in the xz-plane at the moment of emission. Then the direction cosines of the 
non-rotating dipolar magnetic field far from the emission point are 

6,o = <5/2, 6yo-0, 6,o = l-6'o/2- (3.4) 

Now let the magnetic field rotate at the angular velocity 17, the rotational 
axis making the angles a and ^ with the magnetic axis and the wave vector, 
respectively. The open field fine tube is befieved to be narrow, so that |q; — ^| < 
(5 <C 1. To the first order in Qt, the direction cosines of the rotating field are 
then given by 

Ox = Oxo + iltsmt, z , 



by^±ntsm£,—--^, (3.5) 

/)2 -1-7)2 

2 
with t = z/c. The sign of the component by is that of the production b^ • (k x Jl); 
here hd is aligned with the magnetic axis at the moment of emission. Using Eq. 
(3.5) and taking into account that the waves are emitted close to the magnetic 
axis, X, y ^ z, one can obtain the components of the vector q = b x (Jl x r): 

flz . ^a — £ 
q-x == T — sm^^^, 
c d 



qy^ — smC-^^ . (3.6) 

c 



The sign of q^ is opposite to the sign of bc( • (k x fi). 



3.2 Wave mode coupling in case -^ ^ 1 

First we assume that — ^ <^ 1. Then taking into account that ujp oi N oi z^^ 
and substituting Eqs.(3.5), (3.6), one can reduce Eq.(2.8) to the form 



du " "^ r^S 



iuUx = i — ^ri{biuax ± &2ay), 



pi = ±t^nb,ax, (3.7) 

where u = Zp/z, with Zp referring to the polarization-hmiting radius, rj = 

2sin^, 6i = — T ^^, 62 = — T-^. Above we used the definition of the 

polarization-hmiting radius (Eqs. (1.1), (1-2)) implying that 

8ujl{zp)zp ^ ^ 

The wave mode coupling caused by the magnetosphere rotation is described by 
the right-hand sides in Eq. (3.7), which are small as — ^. It is easy to find the 

solutions of Eq. (3.7) to the leading order in — ^. If, say, the ordinary mode is 
emitted, that is at the emission origin Ux = Cx, ay ~ 0, the wave polarization 
can be found to vary along the trajectory as follows: 

ax = Cj;exp(m^/2), 

z /■" 

ay ^ ±i^T]b2Cx / exp(m'V2)dM'. (3.8) 

Setting z ^ 00 (u ^ 0) yields the limiting polarization of the emitted ordinary 
wave: 

ax ^ ^x^ 



.^^^^Cx^^rjb,^ 



exp(i7r/4). (3.9) 



In case of extraordinary wave (at emission ax — 0, ay — Cy), the evolution of 
polarization can be treated similarly. The limiting polarization is then given by 

Z I IT 

ax = TiCy—^-qb2J-cyiY>{-iT^/'^), 

ay^Cy, (3.10) 



So the wave mode coupling in the polarization-hniiting region causes the elhp- 
tical polarization of the escaping normal waves. The contribution of circular 
polarization is characterized by the normalized Stokes parameter V defined as 

ao,ai + ay a* 
Involving Eqs. (3.9), (3.10) one can obtain 

|y|=4^sine^. (3.12) 

rLO 

It is easy to see that the handedness of the circular polarization resulting from 
the limiting polarization effect remains constant throughout the pulse. Note 
that, to the first order in — ^, the polarization profile appears to be symmetri- 
cal, with the peak at the centre. This agrees qualitatively with the observational 
data for a number of pulsars (Han et al. 1998). According to Eq. (3.12), the de- 
gree of circular polarization is ^ — ^ . Although this quantity was regarded as a 

2: 

small parameter, it proved to be not very small at pulsar conditions, — ^ > 0.1 

(Cheng & Ruderman 1979). Thus, the linearly polarized waves can acquire 
considerable circular polarization because of the wave mode coupling in the 
polarization-limiting region. In particular, the highest observed circular polar- 
ization, ~ 60% for PSR 1702-19 (Biggs et.al. 1988), is likely to be explained by 
this effect. 

3.3 The case ^ > 1 

Now wc turn to the wave mode coupling in the opposite limit, — ^ ^ 1. The lat- 
ter inequality corresponds to short-period pulsars (Barnard 1986). Substituting 
Eqs. (3.5) and (3.6) into Eq. (2.8) we have, to the first order in '^-^'■ 

dai i i tlS ^u 3i tlS , 1/4 

-1 7^1 = Tt^ ^-zb2W ' 02 - — ^-biw ' oi, 

aw 4 loZpSm^ 8 ZpSm^ 

=F — ^—b2W ' oi, (3.13) 



dw 16zpSin^ 



with 



fli = feifla, - &2ay, a2 ^b2ax + biay, w^{zp/z). (3.14) 

Here the polarization-limiting radius, Zp, is given by the relation: 

ucy^Zp sin^ ^ 



Solving Eq. (3.13) to the leading order in ^-^, one can trace the evolution of 
the polarization of the ordinary wave (initially ai — Ci, a2 — 0)' 

ai = Ci exp(zw/4), 

rLSb2 



a2^TCi •"""" 4^1/4 exp(iu;/4) - / w'--^'^ ejip(iw' /4)dw' . (3.15) 
16zpSin4 I J J 

Setting w — > we find the limiting polarization: 

ai =Ci, 

_r< rLSb2 r(i/4) 

a2 = TCi— ^^TTTTTTTiexp i7r/8 . 3.16 

iDZpSm^ (1/4)^/* 

Correspondingly, the limiting polarization of the extraordinary wave (initially 
polarized as ai — 0, a2 — C'2) can be found to be 

, ^ rL5b2 r(l/4) . 

iD^pSm^ (1/4)^/^ 

a2^C2. (3.17) 

With Eqs. (3.16) and (3.17), the normalized Stokes parameter becomes 

1^1^0.2 "^^'"7^' . (3.18) 

Zp sm ^ 

Now the circular polarization is constant along the pulse, with \V\ < 10%. Such 
values of V are observed for a number of pulsars. Note that in this case the 
handedness of circular polarization also remains constant throughout the pulse. 



4 Conclusion 

We investigated the evolution of normal wave polarization along the trajec- 
tory in the ultrarelativistic highly magnetized electron-positron plasma filling 
the magnetosphere of a pulsar. The plasma allows two normal waves which 
can ultimately escape from the magnetosphere. In the emission region, they are 
linearly polarized, the radiation emitted being also linearly polarized or contain- 
ing an incoherent mixture of radiation with the two linear polarizations. While 
propagating inside the magnetosphere, the waves are influenced by the plasma. 
So the emergent polarization should be set up by the propagation effects. 

The polarization behaviour was described proceeding from Maxwell's equa- 
tions together with the equation of particle motion in the rotating magneto- 
sphere. For simplicity we examined the wave mode coupling introduced by the 

z z 

rotation of a dipolar magnetic field in the limits — ^ <C 1 and — ^ ^ 1, where 

10 



Zp is the polarization-limiting radius, r^ the light cylinder radius, 5 the beam 
width. As a result of this effect, the emergent radiation should be elliptically 
polarized, with the sense of the circular component being the same throughout 
the pulse. Given that r^d <^ Zp < r^, the degree of circular polarization is 
found to be of a few per cents and to remain constant along the pulse. In case 
Zp < T]^5 the wave mode coupling can cause considerably higher circular polar- 
ization. In this limit, the polarization profile is symmetrical and peaks at the 
center of the pulse. Such behaviour is compatible with the observed V-profiles 
in many pulsars. So the observations testify to the case Zp < tlS implying that 
the polarization-limiting radius should lie not too far from the emission region. 
This is also supported by the absence of essential decrease in the position angle 
swing as a result of the limiting polarization effect (Lyne & Manchester 1988). 
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